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peak powers, and show that other users do not transmit or transmit at their corresponding peak power. 
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I. Introduction 

Cognitive radio is a promising technology for improving spectrum utilization in wireless 
communications systems [HI. A secondary user (SU) in a cognitive radio network is allowed to 
access the licensed spectrum allocated to a primary user (PU) if the spectrum is not utilized by the 
PU. Such a spectrum sharing strategy, which is referred to as spectrum overlay or opportunistic 
spectrum access (OS A) [|2l, requires correct detection of spectrum opportunities by the SU. 
Existing works on spectrum overlay have mainly studied spectrum sensing and access policies 
at the medium access control (MAC) layer [j3l- (H. An alternative strategy, which is known as 
spectrum underlay [|7l- [[T2|. enables the PU and the SU to transmit simultaneously, provided 
that the received interference power by the PU is below a prescribed threshold level. A number 
of works have recently studied information theoretic limits for resource allocation in the context 
of spectrum underlay. 

In [fT3|. the optimal power allocation which aims at maximizing the ergodic capacity achieved 
by an SU is derived for various channel fading models subject to the peak interference power 
(PIP) constraint or average interference power (AIP) constraint imposed by a PU. In [fI5l, 
the authors derive the optimal power allocation for the ergodic capacity, outage capacity, and 
minimum-rate capacity of an SU under both the PIP and AIP constraints from a PU. The ergodic 
capacity, delay-limited capacity, and outage capacity of an SU is studied in [[TSl under different 
combinations of the peak transmit power (PTP) constraint or average transmit power (ATP) 
constraint at the SU and the PIP constraint or AIP constraint from a PU. However, all the 
papers mentioned above consider the setup of a single SU. The most recent work fT6l| studies 
a cognitive radio network of multiple SUs under multiple access channel and broadcast channel 
models, where the optimal power allocation is derived to achieve the maximum sum ergodic 
capacity of the SUs subject to various mixed transmit and interference power constraints. The 
optimality conditions for the dynamic time division multiple access scheme are also derived. 

In this paper, we focus on a cognitive radio network where multiple SUs share the licensed 
spectrum of a PU using the frequency division multiple access (FDMA) scheme. The sum ergodic 
capacity of the SUs, which is a relevant network performance metric for delay-tolerant traffics, is 
studied. Besides the transmit power constraints at the SUs and the interference power constraint 
imposed by the PU, which are also considered in lfT3l - lfT6l . we also take into account the total 
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bandwidth constraint of the shared spectrum. Such study is motivated by the fact observed for 
a number of different applications that joint bandwidth and power allocation can significantly 
improve the performance of systems with limited both individual (power) and public (bandwidth) 
resources [[T7l - [|23l . Thus, in this paper, instead of conventional fixed and equal bandwidth 
allocation used in FDMA, we investigate dynamic and unequal bandwidth allocation, where 
the bandwidth allocation varies for different SUs at different channel fading states. Moreover, 
different from the existing works lfT3l - lfT6l . all combinations of the transmit power constraints 
and the interference power constraints are considered, including both PTP and ATP constraints 
combined with both PIP and AIP constraints. 

We first derive the optimal bandwidth allocation for any given power allocation in any channel 
fading state, which results in equivalent problems that only involve power allocation. Using 
the convexity of the resultant power allocation problems, we apply dual decomposition which 
transforms these problems into equivalent dual problems, where each dual function involves a 
power allocation subproblem associated with a specific channel fading state. The dual problems 
can be solved using standard subgradient algorithms. For the power allocation subproblem 
under all combinations of the power constraints, we derive the structures of the optimal power 
allocations. These structures indicate the possible numbers of users that transmit at nonzero 
power but below their corresponding peak powers, and show that other users do not transmit or 
transmit at their corresponding peak power. Based on these structures, we develop algorithms 
for finding the optimal power allocations in each channel fading state. 

The rest of the paper is organized as follows. Section |ll] summarizes the system model and 
formulates corresponding sum ergodic capacity maximization problems. Section Hill derives the 
optimal bandwidth allocation for the problems formulated in Section HI] subject to the bandwidth 
constraint. Section |IV] obtains the optimal power allocations from the resultant problems in 
Section Ull] under all combinations of the transmit power constraints and interference power 
constraints. Numerical results for the maximum sum ergodic capacity under different combina- 
tions of the power constraints and the bandwidth constraint are shown in Section |Vl Section |VI] 
concludes the paper. 
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II. System Model 

Consider a cognitive radio network of SUs and one PU. The PU occupies a spectrum of 
bandwidth W for its transmission, while the same spectrum is shared by the SUs. The spectrum 
is assumed to be divided into distinct and nonoverlapping flat fading channels with different 
bandwidth, so that the SUs share the spectrum through FDMA to avoid interferences with each 
other. The channel power gains between the ith SU transmitter (SU-Tx) and the ith SU receiver 
(SU-Rx) and between the ith SU-Tx and the PU receiver (PU-Rx) are denoted by hi and 
respectively. The channel power gains, i.e., g = [gi g2 ■ ■ ■ Qn] and h = [hi h2 ■ ■ ■ h^], are 
assumed to be drawn from an ergodic and stationary vector random process. We further assume 
that full channel state information (CSI), i.e., the joint probability density function (PDF) of the 
channel power gains and the instantaneous channel power gains, are known at the SUsQ The 
noise at each SU-Rx plus the interference from the PU transmitter (PU-Tx) is assumed to be 
additive white Gaussian noise (AWGN) with unit power spectral density (PSD). 

We denote the transmit power of the ith SU-Tx and the channel bandwidth allocated to the 
ith SU-Tx as pi(g,h) and Wi{h,g), respectively, for the instantaneous channel power gains g 
and h. Then the total bandwidth constraint can be expressed as 

N 

Y,Mh,g)<W, Wh,g. (1) 

i=l 

The PTP constraints are given by 

p,{h,g)<Pf\ Wz,h,g (2) 

where Pf denotes the maximum peak transmit power of the ith SU-Tx. The PIP constraint is 
given by 

N 

J29iPi{h,g)<Q^\ ^h,g (3) 

1=1 

where Q^^ denotes the maximum peak interference power allowed at the PU-Rx. The ATP 
constraints are given by 

E{p,ih,g)}<Pr, yt (4) 

'Note that the full CSI assumption is typical in the context of cognitive radio and is also made in other works such as 1131 - 
1161 . Indeed, under this assumption we aim at investigating the information-theoretic limits on the sum ergodic capacity. 



June 29, 2010 



DRAFT 



5 



where the expectation is taken over h and g, and Pf-"" denotes the maximum average transmit 
power of the ith SU-Tx. The AIP constraint is given by 



N 



E\Y.gMh,9)\ <Q'''' (5) 



i=l 



where Q"" denotes the maximum average interference power allowed at the PU-Rx. 

The objective is to maximize the sum ergodic capacity of the SUs, which can be written as 

E I y w.{h, g) log fl + ^ 1 (6) 

{w,(h,g)Mh,9)}&r V Wi{h,g) J j 

where is a feasible set specified by the bandwidth constraints ([T]) and a particular combination 
of the transmit power constraints I©, dH)} and the interference power constraints {dS]), dH)}. Note 
that the constraints on nonnegativity of the bandwidth and power allocations, i.e., Wi{h,g) > 
and Pi{h, g) > 0, Vz, h, g, are natural and, thus, omitted through out the paper for brevity. 

It can be shown that the objective function of the problem ^ is concave with respect to 
{wi{h, g),pi{h, g)}, Wi, h, g. It can also be seen that the bandwidth and power constraints ([T])- 
([5]) are linear and, thus, convex. Therefore, the sum ergodic capacity maximization problem Q 
under different combinations of the constraints ©-© is a convex optimization problem. 

III. Optimal Bandwidth Allocation 

Given that the power allocation Pi{h, g), h, g, is fixed, the maximum sum ergodic capacity 
can be expressed as E{/o(^,sr)}, where fo{h,g) is given by 

N 

fo{h, g) = inax V d {wi{h, g)) (7a) 

{wi{h,g)} ^ 
1=1 

N 

s.t. Y,w,{h,g)<W (7b) 



(7c) 



where Gi{wi{h, g)) = Wi{h, g) log (1 + hiPi{h, g)/wi{h, g)) is an increasing and concave func- 
tion of Wi{h,g). The problem (r7al)- (|7bl) is similar to the classical water-filling power allocation 
problem. Thus, the optimal solution of the problem (TTab -dTbl). denoted by {w-(/i,gr)}, must 
satisfy 

dGi{wi{h,g)) _ dGj{wj{h,g)) ^i^j (8) 



dwi{h,g) 



dGj{wj{h,g)) 



«,.(M)=<(M) dw,ih,g) 



5 

Wj{h,g)=w'Jh,g) 
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Since we have 

dGi{wi{h,g)) 



, ^ hiPi{h,g)\ h,pi{h,g) f hiPi{h,g) 

log 1 + — 7777^ --777—^-^-7777:7 = ^ 



w'Ah,g) = w J;'^:'"''j'J ^ . (11) 



dwi{h,g) ^^^h,ghn'^ih,g) V w'i{h,g) J w'^ih, g) + hiPi{h, g) \w'i{h,g) 

(9) 

where Y{x) = log(l + x) — x/{l + x) is a monotonically increasing function, we can obtain 
from dSl) that 

hip,{h,g) _ hjpj{h,g) , 
w[{h,g) w'j{h,g) 

It follows from (iTbl) that at optimality we have '^^=iw'i{h, g) = W. Furthermore, using (fTOl) . 
we can obtain that 

T.f=ihiPt{h,g) 

Substituting the optimal Wi{h,g) given by (fTTI) into dll), we can equivalently rewrite ([61) as 

max E < VV log 1 + ' > (12) 

{Mh,9)}eT' y y ^ J) 

where J-'' is a feasible set specified only by a particular combination of the power constraints 
{dll), ([3]), dH), ©}. Therefore, the optimal power allocation obtained from the problem © and 
denoted by {p*{h,g)}, can also be obtained by solving the equivalent problem (fT2l) . Then the 
optimal bandwidth allocation obtained from the problem ^ and denoted by {w*{h, g)}, can be 
found as 

M-w^fh^ (13) 

IV. Optimal Power Allocation 

In this section, we study the optimal power allocation obtained from the problem (fT2l) with 
J-'' specified by different combinations of the power constraints. 

A. Peak transmit power with peak interference power constraints 

Consider F' = {the constraints Q and ([3])}. Then the optimal value of the problem (fT2l) can 
be expressed as E{fi{h, g)}, where fi{h,g) is given by 

fi{h,g) 



max ly log ( 1 + 2^1^ 1 (14a) 

{Mh,9)} \ W 



s.t p,{h,g)<P[\ Wz (14b) 



N 



Y.9^P^{h,9)<Q''^. (14c) 



1=1 
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For brevity, we drop the dependence on h and g that specifies instantaneous channel power gains. 
Also let {p*} denote the optimal solution of the problem (I14al) - (ll4cl) . Introducing = giPi, the 
problem (|14a| )- (|14c| ) can be equivalently rewritten as 

maxy^—Qi (15a) 



s.t. q, < g,Pf\ V I (15b) 

N 

J2q^<Q''''. (15c) 

i=l 

Let {q*} denote the optimal solution of the problem (I15al) - (ll5cl) and (si,S2,--- ,sn) denote 
a permutation of the SU indexes such that hg^/gsi > hg^/ > ■ ■ ■ > hgf^/gsi^- It is assumed 
that hi/gi ^ hj/gj, Vi ^ j, since hi, gi, hj, and gj are drawn from a continuous-valued random 
process. Then the following lemma is in order. 

Lemma 1: There exists k, 1 < k < N, such that q* = OgPf^, Wi, 1 < i < k - 1, < q* < 
gs.P^K and q* = 0, Vi /c + 1 < i < iV. 

Proof: Let g*^ > for some j and let / < j for some /. First we prove that g* = gsiPsj' by 
contradiction. If g* < gg^Pf^, then we can always find Ag > and define a feasible solution 
{g^J of the problem ([I5al)-([l5cl) q'^^ = q*^ - Ag, g^ = g^ + Ag, g^^ = q*^,\/i, i^j^i^l such 
that the objective function in (I15al) achieves larger value for {g^.} than for the optimal solution 
{g*}, since we have 

i:^i,-i:^c-(^-^)A,>o. (16) 

\gsi ysj / 

Therefore, it contradicts the fact that {g^.} is the optimal solution of the problem (I15al) - (ll5cl) . 

Let ql. < gsjPg^ for some j and let / > j for some /. Using the result obtained above, it can 
be proved also by contradiction that g*^ = 0. This completes the proof. □ 

Lemma 1 shows that for the optimal power allocation under the constraints Q and ([3]), there 
exists at most one user that transmits at nonzero power and below its peak power, while any 
other user either does not transmit or transmits at its peak power. 

Note that either the constraints (|15bl) or the constraint (I15c|) must be active at optimality. 
Using the structure of {g*} given in Lemma 1, k can be found by Algorithm 1. 
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Algorithm 1 Algorithm for finding k in Lemma 1 
Initialize: k = 1 

while Yli=i 9s.Pt <QP'' and k < N -1 do 

k = k + l 
end while 
Output: k 



Since p*^ = q*Jgs,, we obtain 



Ps 



Pf, l<i<k-l 



min{P^;^ {Q^^ - Eti 9s.Pt)l9s.}. ^ = k (17) 
0, k+l<i<N. 

Note that for brevity, we say in this paper that Yl^=i Xi = if n = with a little abuse of 
notation. 

B. Average transmit power with average interference power constraints 

Consider = {the constraints dH) and ©I. Then the dual function of the problem (fT2l) can 
be written as 

N 

/2({A,}, ^) ^ E 9)} + Y, \Pr + f^Q"'' (18) 

i=l 

where {Aj|l < i < N} and /i are the nonnegative dual variables associated with the correspond- 
ing constraints in dH) and ^ and f^ih^g) is given by 

m,g) 4 max log ( i 9) \ _f2j^p^ih,g) (19) 

with 7j = Aj + fiQi. Let {p*} denote the optimal solution of the problem (fT9l) . where we drop 
the dependence on h and g for brevity. Also let F{{pi}) denote the objective function in (fT9l) . 
If p* > for some i, the following must hold 



dF{{p,}) 



dpi 

Then the following lemma is of interest. 

Lemma 2: If hi < •ji for some i, then p* = 0, 



- 7^ = 0. (20) 
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Proof: If p* = 0, Vj, then = 0. If p* 7^ for some j, it can be seen that (|20l) can not be 
satisfied since hi < 7^. Thus, p* = 0. □ 
If p* = for some i, the following must hold 



= ^ 7i<0. (21) 



dpi 

Then the next lemma is in order. 

Lemma 3: p* = 0, Vi if and only if hi < 7^, Vi. 

Proof: It can be seen from Lemma 2 that if hi < 7j, Vi, then p* = 0, Vi Moreover, it can 
be seen from dlD that if p* = 0, Vi, then hi < □ 

Let (si, S2, ■ ■ ■ , Sat) denote a permutation of the SU indexes such that hs-^/'jsi > hs2hs2 > 
■ ■ ■ > hgj^/'jsN- Then we can also prove the following lemma. 

Lemma 4: There exists at most one k such that pi > 0. Moreover, k = si. 

Proof: We prove it by contradiction. It can be seen from (|20l) that if p* > and P ^ > for 
some i ^ j, the following must hold 

hi hi 

— = ^. (22) 

Since hi, 7^, hj, and 7-, are independent constants given in the problem (fT9l) . ((22l) can not be 
satisfied. Let p^ > and p* = 0, Vi, i ^ k. Then it follows from dJO]) and ([21]) that the following 
must hold 

- > V ^ ^ /t. (23) 

Ik It 

Therefore, we must have k = S\. □ 
Lemma 4 shows that for the optimal power allocation under the constraints dH) and ©, there 
exists at most one user that transmits at nonzero power, while any other user does not transmit. 
Case 1: Consider the case when hi < 7j, Vi. It follows from Lemma 3 that p* = 0, Vi. 
Case 2: Consider the case when hi < 7^ does not hold for some i. Using Lemma 4, let 
pI > and p* = 0, Vi, i ^ k. Substituting {p*} into dSO]), we have p*^ = W{l/^s^ - l/hsj. 
Therefore, we obtain 

. ' W{l/{X,,+i,g,,)-l/hJ, z = l 
Ps. = < (24) 
0, 2<i< N. 



June 29, 2010 



DRAFT 



10 



C. Peak transmit power with average interference power constraints 

Consider F' = {the constraints Q and ©j. Then the dual function of the problem (fT2l) can 
be written as 

/3(/i)^E{/^(/i,5)} + /.g- (25) 

where /x is the nonnegative dual variable associated with the constraint Q, and f^ih, g) is given 
by 

g) ^ ^ inax W log (l + SkMiM^ _ ^ (26a) 

{piCi.s)} y w I ^—^ 

s.tp^{h,g)<Pt, Vz. (26b) 

Let {p*} denote the optimal solution of the problem (|26al) - (|26b|) after dropping the dependence 
on h and g for brevity. The following cases are of interest. 

Case 1: Consider the case when hi < figi. Mi. Since the problem (|26a| )- (|26b|) without the 
constraints (|26bl) has the same form as the problem (fT9l ). and pi = 0, Vi, satisfies the constraint 
(I26bl) . it can be seen from Lemma 3 that p* = 0, Vz. 

Case 2: Consider the case when hi < ^Qi does not hold for some i. The problem (I26al) - (l26bl) 
is equivalent to 

max TV log f 1 + SkMM ) q. (27a) 

s.t. qi < fig.Pt, V t (27b) 

where qi = yUfifiPj. Let {q*} denote the optimal solution of the problem (|27al )- (|27bl ) and 
(si, S2, ■ ■ ■ , stv) denote a permutation of the SU indexes such that hs^/ jJ-Qs^ > hs^/ngs2 > ■ ■ > 
hg^l jjigsf^- Then the following lemma is in order. 

Lemma 5: There exists k, 1 < k < N, such that q*. = gsiPfl', Vi, 1 < « < A; — 1, < g*^ < 
gs^Pf, and q* = 0, Mi, k + l< i < N. 

Proof: Consider the following intermediate problem 

max y — -qi (28a) 
M ~l l^9i 

s.t. q, < fxg,P^\ V t (28b) 

N 

Y,1^ = Q (28c) 

i=l 

June 29, 2010 DRAFT 



11 



where Q = J2iLi1i ^^'^ unknown since {q*} is unknown. Let {q^} denote the optimal 
solution of the problem (|28aM28cl) . If {q-} 7^ {q*}, we have Y^iLi hiq'i/ fJ^Qi > Y^iLi ^il* / f^9i 
since {q*} is a feasible solution of the problem (|28a|) - (|28cl) . Then we have 

F m) - F ({g*}) = log [l + SiMM j _ H^iog + ELMM j > (29) 

where denotes the objective function in the problem (I27al) - (l27bl) . Since {q[} is a feasible 

solution of the problem (I27al) - (l27bl) . it contradicts the fact that {q*} is the optimal solution of 
the problem (I27al) - (l27bl) . Therefore, it must be true that {q[} = {q*}. 

It can be seen from the constraints (|27bl) that Xlili — Sili Qi = Q ^ Sili f^9iFi''- Then 
the problem (I28al) - (|28cl) is equivalent to the following problem 

^ h- 

max y — —qi (30a) 

l^gi 

s.t. q, < fxg,P[\ V I (30b) 

N 

J2q^<Q (30c) 

1=1 

since the constraint (I30cl) is active at optimality. Therefore, the problem (I27al) - (l27bl) is equivalent 
to the problem (|30a| )- (|30cl ). Since the problem (|30a| )- (|30c| ) is similar to the problem (|15a| )- (|15c| ) 
in Section HV-Al we conclude that {q*} has the same structure as that given in Lemma 1. □ 

The result of Lemma 5 is similar to that of Lemma 1 . Specifically, it shows that for the optimal 
power allocation under the constraints ^ and ([5]), there exists at most one user that transmits 
at nonzero power and below its peak power, while any other user either does not transmit or 
transmits at its peak power. 

Using Lemma 5, let q*. = ^.Qs^Pf^, Wi, 1 < i < k - 1, < q*^ < /igf^.P^'', and q*. = 0, Wi, 
k + 1 < i < N . Then we only need to find k and q*^ to determine {q*}. 

Consider the case when < g*^ < jigg^Pg^, 1 < k < N. Then the following must be true 




-1 = (31) 



W 




(32) 
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Substituting {q*J into dlB, we obtain = W{1 - figsjhs,) - l^9s, EiJ KP^/hs,- Since 
q* must satisfy < g* < ^Qs^P"^^, it must be true that 



k-l / , X k 



Y: K.Pt < (-^ - l) < E (33) 



j=l ^ ' i=l 

ok 



Consider the case when g*^ = ^gs^P^^, 1 < k < N — 1. Then the following must hold 

dH{qs, 



dqs 

and 



dq 



K+i/l^9. 



Sfe + 1 



1 < 0. (35) 



(36) 



Substituting {q*} into dH) and ([35]), we obtain 

W ( _ 1^ < y hs-Pt <w(^-lY l<k<N-l. 

\^^9s,+, J ^ ' ' \f^9s, J 

If q*^ = jjigs^P"^^, k = N, then only (|34l) must be true and it follows that 

y /i,,Pf < W - iV A; = iV. (37) 

Lemma 6: T/zere ejcz5?5 on/y one set of values for {q*} that satisfies only one of the necessary 
conditions dSB, (EH) or (l35l) . 

Proof: It is equivalent to prove that there exists only one k that satisfies only one of (|33l) . 
or dlT]). Let = ELi - ^(^^,//^^^, - 1) for brevity. Then it must be 

true that Lq < Li < ■ ■ ■ < Ln, Mi > M2 > ■ ■ ■ > Mn and Lq < Mi. It can be seen that if dST]) 
holds, i.e., if Li < Mi, ^i, 1 <i< N, then ^ and ^ do not hold. 

If (l37l) does not hold, then these exists such / that Li < Mi, yi, 1 < i < I ~ 1 and Li > Mi, 
yi, 1 < i < N. The following two cases should be considered, (i) If Li_i < Mi < Li, (l33l) 
holds for k = I. Since Li < Mi, Vi, 1 < z < / — 1, (l33l) does not hold for k < I as well. Since 
Mi < Ml < Li < Vi, / + 1 < i, (I33l) does not hold for /c > /. Since Li < Li+i < M^+i, 

yi, 1 < i < I - 2, dlH) does not hold for A; < / - L Since < Mi, dlH) does not hold also 
for A; = / — L Moreover, since Mj < Li, Mi, I < i, (l36l) does not hold for A; > / — L Therefore, 
only (133]) holds for only k = I. (ii) If Mi < Li^i < Mi^i, ^ holds for A; = Z - L Similar to 
the case (i), it can be proved that only (l36l) holds for only k = I — 1. □ 
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Algorithm 2 Algorithm for finding k in Lemma 5 
Initialize: A; = 0, c = 

while c = do 

k = k + l 

if Eti KPt < WihJ^^g,, - 1) < Ell hs^Pt then 

c = 1 
end if 

if {Wihs.Jfigs,^, - 1) < Ell hs.Pl' < WihsJ^^g,, - 1) and A: < AT - 1} or 

{Ell hs.Pt < Wihsjf^gs, - 1) and k = N} then 
c = 2 

end if 
end while 
Output: k, c 



Using Lemma 6, Algorithm 2 is developed to find the unique k in Lemma 5. Note that k 
satisfies (l33l) and (l36l) or (l37l) if the output of Algorithm [21 is c = 1 and c = 2, respectively. 
Since p*. = q*J figsi, when c = 1, we obtain 



Ps 



1 < z < A; - 1 

1^(1//.^?,, - lA.J - El>s,P.f As,, ^ = A; , l<z<Ar (38) 

0, A; + 1 < i < 



and when c = 2, we obtain 



Ps 



Pf.^ l<i< k 



0, 



A; + 1 < z < A^ 



, 1 < i < A^. 



(39) 



D. Average transmit power with peak interference power constraints 

Consider = {the constraints © and Then the dual function of the problem (fT2l) can 
be written as 

N 

U{X,})^E{j2h,g)} + Y,\Pr (40) 



i=l 
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where {Aj|l < i < N} are the nonnegative dual variables associated with the corresponding 
constraints dH) and f^{h,g) is given by 

9) ^ ^ max W log (l + E^^li^^ _ ^ A.p.(^, g) (41a) 



N 



s.t. 5^(?.p.(^^)<Q^^ (41b) 



i=l 



Let {p*} denote the optimal solution of the problem (|41al) - (l41bl) where the dependence on h 
and g is dropped for brevity. The following three cases are of interest. 

Case 1: Consider the case when hi < Aj, Vi. Similar to Case 1 in Section |IV-C[ it can be 
seen from Lemma 3 that p* = 0, Vi. 

Case 2: Consider the case when hi < \i does not hold for some i and the constraint (I41bl) 
is inactive at optimality. Let (si, S2, ■ ■ ■ ,sn) denote a permutation of the SU indexes such that 
hg-^/Xg^ > hg^/ > ■ ■ ■ > hgj^/Xgj^. Since the problem (I41al) - (l41bl) without the constraint (I41bl) 
has the same form as the problem (fT9l) . it can be seen from (l24l) that p*^ = W{l/\s-^ — l/ZisJ 
and p*^ = 0,\fi, 2 < i < N, if it satisfies the constraint (I41bl) . i.e., YliLi 9s,P*s, = 9siW{l/Xs^ - 
1/hsJ < QP'. 

Case 3: Consider the case when hi < Aj does not hold for some i and the constraint (|41bl) is 
active at optimality, i.e., gs^W{l/ Xg^ — ^/hgj > Q^'^. The dual function of the problem (I41al) - 
(I41bl) can be written as /4 (/x) = fl' + H-Q^^, where ^ is the nonnegative dual variable associated 
with the constraint (|41bl) . and f'l' is given by 

r: 4 maxl^log 1 + ^'^=\p^' - 5^ A.P, -^^Y.g.p.. (42) 

Let /i* denote the optimal dual variable. Also let F{{pi}) denote the objective function in the 
problem (|42|) . If p* > for some i, the following must hold 



dF{{pi}) 



dpi 



a/V— T777 -^i- ^^*9^ = 0. (43) 



If p* = for some i, the following must hold 

dFiip,}) 



dpi 



^\ - A. - ^l*g^ < 0. (44) 



Note that since the problem (I41a|) - (|41bl) is convex, the necessary conditions (l43l) and (l44l) for 
the optimal solution {p*} are also sufficient conditions. 
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Lemma 7: There exists at most two j k such that p* > and pi > 0. 

Proof: We prove it by contradiction. It can be seen from (|43l) that if p* > 0, p* > 0, and 

pI > for some i ^ j, j ^ k, i ^ k, the following must hold 

hi _ hj _ hk ^^^^ 



Xi + fi*gi Xj + fi*gj Xk + fi*gk 
Since hi, Xi, gi, hj, Xj, gj, hk, Xk, and gk are independent constants given in the problem 

(I41al) - (l41bl) . and only fi* is a variable, (l45l) can not be satisfied. □ 
Lemma 7 shows that for the optimal power allocation under the constraints ([3]) and dH), there 

exists at most two users that transmit at nonzero power, while any other user does not transmit. 
Then Case 3 can be further divided into the following two subcases. 

Case 3.1: Consider the subcase when pi > and p* = 0,\/i ^ k. Since the constraint (I41bl) is 
active at optimality, i.e., J^iLi 9iPi = 9kPl = Q^^^ we obtain that pi = Q^^ /gk- Then substituting 
{p*} into (|43l) we have 

^* ^ Qk/hk + Qr'^/W ' Vk ^"^^^ 
Note that /i* given in (|46l) must satisfy /i* > 0. Substituting {p*} into (|44l) . we can see that /i* 

given in (|46l) also must satisfy 

1 + hkQP^/gkW gi 

Then Algorithm 3 can be used to find k. Note that {p*} does not exist in Case 3.1 if the output 
Algorithm 3 Algorithm for finding k in Case 3.1 



k = argmax{i} W\og (l + 



1 Afe 



if /i* < max{,^,| i+hmI'),,w - | or /.* < then 

= 
end if 
Output: k 



of Algorithm [3] is k = 0. 

Case 3.2: Consider the subcase when p* > 0, pi > 0, j ^ k and p* = 0, Vi, i j,i k. It 
follows from (gS]) that 

hj _ hk 



Xj + fi*gj Xk + II* gk ' 



(48) 
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Therefore, we obtain that 

^* = - Xk/hk ^ ^^^^ 

Qk/hk - Qj/hj ■ 

Note that n* given in (l49l) must satisfy /U* > 0. Using (l43]) and the fact that the constraint (I41bl) 
is active at optimality, we have 

+ hkPl = Wh,/{\, + - W ^^^^ 
9jP*j + 9kPl = Q^^- 



Solving the system of equation (1501) . we obtain 

. ^ Q^'/gk - a/hk , ^ a//i, - QP'^/gj ^^^^ 
9j/9k - hj/hk ' /ife//ij - Qk/gj 

where a = Whj/{Xj + ii*gj) — W. Note that p* and given in (|5T]) must satisfy p* > and 

> 0. Substituting {p*} and /i* into (|44)) . we can see that j and A; must satisfy 

^4^^^ > VM^J,^^^. (52) 

5'fe/"'fc - gj/hj gi/hi - Qj/hj 

Then Algorithm 4 can be used to find j and A;. Note that {p*} does not exist if the output of 
Algorithm m is j = and A; = 0. 

E. Combinations of more than two power constraints 

Consider J"' = {the constraints ©, and ©} or J"' = {the constraints ©, (g]), and ©}. 
It can be shown that the corresponding dual functions of the problem (fT2l) under these two 
combinations of the power constraints have the same form as those in Subsections IIV-CI and 
IIV-DI respectively. Therefore, optimal solutions can be found similarly therein and, thus, are 
omitted here. 

Consider J"' = {the constraints ©, ©, and ox 7' = {the constraints ©, ©, and ©} 
or J^' = {the constraints @, ([4]), and ©}. It can be shown that the corresponding dual 
functions of the problem (fT2l) under the first two combinations of the power constraints have 
the same form as that under the third combination. Therefore, we only focus on the case 
J-'' = {the constraints©, ©, ©, and ©}. Then the dual function of the problem (fT2)) can be 
written as 

N 

/5({A,}, /i) ^ E {/^(/i, g)} + Y, \Pr + f^Q'''' (53) 



i=l 
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Algorithm 4 Algorithm for finding j and k in Case 3.2 
Initialize: 1=0 

for J = !,■■■ ,N -1 do 

for = j + 1, ■ ■ ■ ,N do 

if /i* > then 

a = Whj/{Xj + fx*g^) - W 

gj/gk-hj/hk ' ^fc H/hj-Qk/gj 
if p* > and > then 

X = XU{(j,A;)} 

= H^log (l + M±;!^) _ x,p* - X,pl 

end if 
end if 
end for 
end for 

(j, k) = argmax{(j,i)gx} Vi^i 

gk/hk-gj/hj i*7^J>«} Qi/hi-gj/hj 

j=0,k = 
end if 

Output: j, A; 



where {Aj|l < i < N} and fi are the nonnegative dual variables associated with the correspond- 
ing constraints in dH) and ^ and f^{h,g) is given by 

l + ^^^^^^P^]-J2^Mh,9)-^^J29^Mh,g) (54a) 



N 



s.t. J29M^^9)<Q'"' (54b) 



2=1 



p,{h,g)<P[\ Vz. (54c) 

Let {p*} denote the optimal solution of the problem (|54a|) - (|54c|) where the dependence on h 
and g is dropped for brevity. The following cases are of interest. 
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Case 1: Consider the case when /ij < Aj + Vi. Similar to Case 1 in Subsections IIV-CI and 
IIV-DI it can be seen from Lemma 3 that p* = 0, Vi. 

Case 2: Consider the case when hi < \i + figi does not hold for some i and the constraint 
(|54b| ) is inactive at optimality. Since the problem (|54a| )- (|54cl ) without the constraint (|54b| ) has 
the same form as the problem (I26al) - (l26bl) . {p*} can be found using Algorithm |2] and (|38l) or 
(l39l) if it satisfies the constraint (I54bl) . 

Case 3: Consider the case when hi < Aj + figi does not hold for some i and the constraint 
(I54bl) is active at optimality. The dual function of the problem (I54al) - (l54cl) can be written as 
= /g" + where (3 is the nonnegative dual variable associated with the constraint 

(|54bl) and f^' is given by 



E«=i hiP 




fi' ^ max log + ^^^^p^ 1 - > : 1^P^ - /3 > : g^P^ (55a) 
s.t. Pi < Pf^ V i. (55b) 

where 7^ = Xi + fxgi. Let /3* denote the optimal dual variable and F{{pi}) stands for the objective 
function in the problem (I55al) . If > p* > for some i, the following must hold 



dF{{p,}) 



dpi 



J"', . 13* 9^ = 0. (56) 



If p* = pf for some i, the following must hold 



dF{{p,}) 



dpi 

Moreover, if p* = for some i, the following must hold 



= aT^ li- P*9i > 0. (57) 



dFiip,}) 



dpi 



Jp^^—777:-l^-^3*9^<0. (58) 



te}=K} 1 + Er=l h^P*/W 

Note that since the problem (I54al) - (l54cl) is convex, the necessary conditions (l56l) . (l57l) and (l58l) 
for the optimal solution {p*} are also sufficient conditions. 

Lemma 8: There exists at most two j and k, j ^ k such that Pj^ > p* > and P^^ > > 0. 

Proof: We prove it by contradiction. It can be seen from dSll) that if Pf > p* > 0, Pf > 

p* > 0, and Pj^'^ > pI > for some i ^ j, j ^ k, i ^ k, the following must be true 



hi _ hj _ hk 



(59) 
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Since h^, gi, hj, 7^, gj, h^, 7a,., and gf, are independent constants given in the problem 
(|54al) - (l54cl) . and only f3* is a variable, (|59| ) can not be satisfied. □ 

Lemma 8 shows that for the optimal power allocation under the constraints ([3]), © and 
©, there exists at most two user that transmit at nonzero power and below their peak power, 
while any other user either does not transmit or transmits at its peak power. 

Then Case 3 can be further divided into the following two subcases. 

Case 3.1: Consider the subcase when Pf^'' > pi > and p* e {Pj''', 0}, Vi ^ k. Let M and 
A/'o denote the sets of SU indexes such that p* = P[^, \fi e Mi and p* = 0, Vi G Mq. Since 
the constraint (I54bl) is active at optimality, i.e., Yl^=i 9iPi = 9kPl + Y^ieMi S'-^i^ ^ Q^''' 
obtain pi = {Qp^ - Y^iaN^ giPt)/gk- Note that pi given here must satisfy Pf > pl>0. Then 
substituting {p*} into (|56l) we obtain 

13* = , ^ , ^. (60) 

1 + [hkiQp'^ - 9^Pt)/gk + h^Pt) /W 9k 

Note that l3* given by ^ must satisfy /S* > 0. Substituting {p*} into dST]) we can see that jS* 
given by (l60l) must satisfy 



r < 7 ^ ^ ^,yteMi. (61) 

1 + [h{QP' - E.eM 9^Pt)/gk + E.eM h^Pt) /w 9^ 

Substituting {p*} into (|58l) . we can see that [5* given in (|60l) also must satisfy 

fi*> 7 ^^^7^ ^, V.GATo. (62) 

1 + [hk{QP' - E.eM 9^Pt)/9k + E.eM ) 

Let ■ ■ denote all the subsets of the set M\{i} where \ denotes the set 

difference operator. Then Algorithm 5 can be used to find k. Mi, and Ao- Note that {p*} does 
not exist if the output of Algorithm [5] is A; = 0. 

Case 3.2: Consider the subcase when Pf > p* > 0, P^'' > > and p* G {Pf'', 0}, 
Vi 7^ j, k. Let Ml and Mq denote the sets of SU indexes such that p* = Pf \/i G Mi and 
p* = 0, Vi G Ao, respectively. It follows from (|56l) that 

^ - (63) 



Therefore, we obtain that 



gk/hk- gj/hj 
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Algorithm 5 Algorithm for finding k. Mi, A/q in Case 3.1 



Initialize: X = 

for A; = 1,2,-- - ,Ar do 
for / = 1,2,--- ,2^-1 do 

M = 5P 

if Pf >pl>0 then 

X = XU {/} 

ri = 
end if 
end for 



W\og (l + ^!^I!^±%ii!!^^ - ^,pi - ^^^^^ 



Vk = max{jgi}ri, t = arg max{jgi} 



^5,* = 5^ 



1 = 
end for 

A; = arg max{j}, f « 

M = SI 

Afo = MWi\{k} 

13* 



hkhk Ik 



l+(?»fc(Qf'^-E,6AAi 9.-Pf' )/9'=+E,:eAAi h,Pf'')/W 9k 

if /3* < or /3* > / , , ^ — 

or /3* < ,^/, ^ , pp.N,,,, - ^, 32 G A/'o then 

A; = 
end if 

Output: k. Ml, A/'o 



Note that (3* given in (l64l) must satisfy (3* > 0. Following (l56l) and the fact that the constraint 
(|54b| ) is active at optimality, we have 

h,P* + hpl = Wh,/{^, + (3*g,) -W- Y..^j^^ KPf (^5^ 
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Solving the system of equation (|65l) . we obtain 

_ a/9k-b/hk ^ ^ b/hj-a/gj ^^^^ 

where a ^ Q^'^ - Y^^^j^, 9rPf and b ^ Wh^jiri, + ii*g,) -W - E.^M ^^^f Note that p* and 

given in ([661) must satisfy Pj'' > p* > and > pl > 0. Substituting {p*} and (3* given 

by dill) into dSH), we obtain 

li/hj--fk/hk ^ lj/hj--fi/hi (g^^ 
gk/hk-gj/hj ~ gi/hi-gj/hj' 

Moreover, substituting {p*} and /3* given by (|64l) into (l58l) . we also obtain 

Ij/hj-Jk/hk lj/hj--ii/h, ^ ■ ^ f,r t.a. 

Let S^^j ,S'fj , ■ ■ ■ , iSfj ■* denote all the subsets of the set Af\{i,j}. Then Algorithm 6 can be 
used to find j, k, A/i, and A/q. Note that {p*} does not exist if the output of Algorithm [6l is 
j = and A; = 0. 

V. Simulation Results 

Consider a cognitive radio network which consists of one PU and four SUs. For simplicity, 
we assume that only Rayleigh fading is present in all links. The variance of the channel power 
gain is set to = 1. We also set W = 1, Pf = 10, Wi, P^" = 10, Wi, = 1, and Q"'' = 1 as 
default values if no other values are specified otherwise. The AWGN with unit PSD is assumed. 
We use 1000 randomly generated channel power gains for h and g in our simulations. The 
results are compared under the following five combinations of the power constraints: the PTP 
with PIP constraints (PTP+PIP), the PTP with AIP constraints (PTP+AIP), the ATP with PIP 
constraints (ATP+PIP), the ATP with AIP constraints (ATP+AIP), the PTP and ATP with PIP 
and AIP constraints (PTP+ATP+PIP+AIP). 

First, we aim at showing by Fig. \T\ that the information-theoretic limit for the sum ergodic 
capacity is indeed significantly higher when bandwidth is allocated optimally as compared to the 
case when it is allocated equally among SUs. In this figure, OBPA stands for optimal bandwidth 
and power allocation, while EBPA stands for equal bandwidth and power allocation. The case 
of PTP+PIP is only depicted in Fig. [H but the conclusion about the superiority of optimal 
bandwidth and power allocation holds true for other combinations of power constraints. Then 
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Algorithm 6 Algorithm for finding j, k. Hi, Afo in Case 3.2 



,pk 
i 



Initialize: X = 

for j = 1, 2, • • • , AT - 1 do 
for ^ = j + 1, • • • , A^ do 
forZ = l,2,-- - ,2^-2 do 

o* _ 'Tj/hj-'Yk/hk 
^ gk/hk-gj/hj 

if /3* > then 

« = Q^' - E.eM ^^^f ' ^ = Wh,/{^, + p*g,) -W- hiPt 

* _ a/gk-b/hk * _ b/hj-a/gj 
9j/gk-hj/hk^ hk/hj-Qk/gj 

if Pf > p* > and Pf > Pfc > then 
X = XU{Z} 

n = wio, (i + M±^!^^^^W^) _ _ _ ^^^^^ ,^p. 
end if 
end if 
end for 

Vj^k = max{jgx} ^ = argmax{jex} n 

c* c(*) 

'^j,k - '^j,k 

end for 
end for 

{j,k) = argmax{(j,/)}^;j,/ 

"/j/hj--/k/hk ^ 7j/hj~fi/hi ■yj/hj-'Yk/hk ^ 1j/hj-"/i/hi »r ^j^^^ 

Qk/hk-gj/hj gi/hi-gj/hj ' ^ gk/hk-gj/hj gi/hi-gj/hj ' 

J = 0, = 
end if 

Output: j. A;, A^o 
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Fig. [2] shows and compares the maximum sum ergodic capacity under PTP+PIP, PTP+AIP and 
PTP+ATP+PIP+AIP constraints versus PP^ where PP'' = Pf, Vi is assumed. It can be seen 
from the figure that the maximum sum ergodic capacity achieved under PTP+AIP is larger than 
that achieved under PTP+PIP for any given P^^. This is due to the fact that the AIP constraint 
is more favorable than the PIP constraint from SUs' perspective, since the former allows for 
more flexibility for SUs to allocate transmit power over different channel fading states. It is also 
observed that the performance under PTP+ATP+PIP+AIP is very close to that under PTP+PIP 
that is because the PTP constraint dominates over the ATP, PIP, and AIP constraints for all 
values of P^^. 

Fig. [3] shows the maximum sum ergodic capacity under ATP+PIP, ATP+AIP and PTP+ATP+ 
PIP+AIP constraints versus P"" where = P^^, Vz is assumed. The maximum achievable 
sum ergodic capacity achieved under ATP+AIP is larger than that achieved under ATP+PIP for 
all values of P'^^ since the PIP constraint is stricter than the AIP constraint. The sum ergodic 
capacity under PTP+ATP+PIP+AIP is much smaller than that under ATP+PIP and ATP+AIP 
due to the fact that the PTP constraint is dominant over other constraints for all values of P'^'" . 

Fig. m shows the maximum sum ergodic capacity under PTP+PIP, ATP+PIP and PTP+ATP+ 
PIP+AIP constraints versus Qp^. It can be seen from the figure that the maximum sum ergodic 
capacity achieved under ATP+PIP is larger than that achieved under PTP+PIP for any given 
Qpk^ This is because the power allocation is more flexible for SUs under the ATP constraint 
than under the PTP constraint. The sum ergodic capacity under PTP+ATP+PIP+AIP saturates 
earlier than that under PTP+PIP and ATP+PIP, because it is restricted by the AIP constraint. 

Fig. [5] shows the maximum sum ergodic capacity under PTP+AIP, ATP+AIP and PTP+ATP+ 
PIP+AIP constraints versus Q"^'. Due to the same reasons as for the results in Fig. |4l the sum 
ergodic capacity achieved under ATP+AIP is larger than that achieved under PTP+AIP. The 
sum ergodic capacity under PTP+ATP+PIP+AIP saturates earlier than that for PTP+AIP and 
ATP+AIP because of the presence of the PIP constraint. 

Finally, Fig. [6] shows the maximum sum ergodic capacity under PTP+PIP, PTP+AIP, ATP+PIP, 
ATP+AIP and PTP+ATP+PIP+AIP versus W . Similar performance comparison results as in the 
previous figures can be observed. One difference is that the sum ergodic capacities do not saturate 
with the increase of W . 
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VI. Conclusion 

A cognitive radio network where multiple SUs share the licensed spectrum of a PU using the 
FDMA scheme has been considered. The maximum achievable sum ergodic capacity of all the 
SUs has been studied subject to the total bandwidth constraint of the licensed spectrum and all 
possible combinations of the peak/average transmit power constraints at the SUs and interference 
power constraint imposed by the PU. Optimal bandwidth allocation has been derived in each 
channel fading state for any given power allocation. Using the structures of the optimal power 
allocations under each combination of the power constraints, algorithms for finding the optimal 
power allocations have been developed too. 
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Fig. 1. Sum ergodic capacity vs P^^ . 
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Fig. 2. Sum ergodic capacity vs P^^ . 
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Fig. 5. Sum ergodic capacity vs Q°". 




Fig. 6. Sum ergodic capacity vs W. 
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